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ABSTRACT 

We present high- and medium-resolution spectroscopic observations of the cataclysmic vari- 
able BF Eri during its low and bright states. The orbital period of this system was found to be 
0.270881(3) days. The secondary star is clearly visible in the spectra through absorption lines 
of the neutral metals Mg I, Fe I and Ca I. Its spectral type was found to be K3±0.5. A radial 
velocity study of the secondary yielded a semi-amplitude of K2 = 182.5 ± 0.9 kms -1 . The 
radial velocity semiamplitude of the white dwarf was found to be K\ = 74 ± 3 kms -1 from 
the motion of the wings of the Ha and H/3 emission lines. From these parameters we have 
obtained that the secondary in BF Eri is an evolved star with a mass of 0.50-0.59 M Q , whose 
size is about 30 per cent larger than a zero-age main-sequence single-star of the same mass. 
We also show that BF Eri contains a massive white dwarf (Mi 1.23Af Q ), allowing us to 
consider the system as a SN la progenitor. BF Eri also shows a high 7-velocity (7 = —94 
kms -1 ) and substantial proper motion. With our estimation of the distance to the system 
(d w 700 ± 200 pc), this corresponds to a space velocity of ^350 kms -1 with respect to the 
dynamical local standard of rest. The cumulative effect of repeated nova eruptions with asym- 
metric envelope ejection might explain the high space velocity of the system. We analyze the 
outburst behaviour of BF Eri and question the current classification of the system as a dwarf 
nova. We propose that BF Eri might be an old nova exhibiting "stunted" outbursts. 

Key words: methods: observational - accretion, accretion discs - binaries: close - stars: 
dwarf novae - stars: individual: BF Eri - novae, cataclysmic variables 



1 INTRODUCTION firmed the classificat ion of BF Eri as a CV of an unknown type 

JChisholm et all 19991) . 

Cataclysmic Variables (CVs) are close interacting binaries that con- „, , ,. , , „ , , , rrrr~ t—\ 

, c The long term light curve was first analysed by Watanabe 

tain a white dwarf (WD) accreting material from a companion, usu- HT^k , j 1 1 1 1111 r.^ ^ • • 

( 1999) who detected two outbursts and concluded that BF En is a 
ally a late main-sequence star. CVs are very active photometrically, ., . , r .» .„ , , , , 

J n J DN with a recurrence period of 40-50 days and the observed map- 

exhibiting variability on time scales from seconds to centuries (see , , „ . . | T , . 1 nr^r 

°l 1 J 1 mtude range between 13.2 and fainter than 14.7 mag. Kato ( 199 

review by Warner 1995). Dwarf novae (DNs) are an important sub- ,rn— - — r-rz — -1 „ , . . , . , 



set of CVs. Observable features that distinguish DN from other 
CVs such as nova-like variables (NLs) are the recurrent outbursts of 



and lKato & Uemural fcOOd) confirmed this classification but men- 
tioned the low amplitude of its outbursts and the difference in the 
outburst activity between seasons as the amplit udes of out bursts 



2-6 mag that occur on timescales of days to years. In this paper we , ,, ,„ , , , \~ IJ,r>„ri 

were remarkably different - 1.0 mag and 1.6 mag. Kato ( 1999) was 
concentrate on BF Endam, a little-studied dwarf nova candidate. „ , . . .... . , 

unable to detect any orbital variability. 
BF Eri was discovered and first classified as a slowly vary- „ . i( , - ... , . , 

I — I 1 — 7-1 . , „ . Prior to the beginning of this work there was no other study 

ing variable star by Hanley & Shaplev (1940). During the fol- JT1 ^, ■ . j r 

~ . — ~~ — of the system in the literaturaj. This motivated us to perform time- 
lowing half a century no observations of the star were reported. , , - r . . , . , .. ... 

J r resolved spectroscopy of BF En in order to study its properties in 

As a consequence, the specific nature of the star remained un- , . ., T ... ... ,, . ,,. , , 

. ' r more detail. In this paper we present and discuss the first high- and 
known. The identification of the Einstein survey source 1ES 0437- .. ,. ., . . , , 

1 ; 11 -1 medium-resolution spectroscopic observations of BF En during its 

046 with BF En (Elvis et al. 1992) and optical spectroscopy by , , , . . of .„, , T „ . pn , 

M ' 1 " ^ . low and bright states in 2005 and 2006. In Section |2| we describe 

19961) eventually led to its correct identification as 

a CV. Later, the ROSAT identification of the Einstein source con- 

1 ISheets et alj l2007t) have recently presented optical medium-resolution 
spectroscopy of BF Eri, obtained primarily to determine the orbital period. 
* E-mail: vitaly.neustroev@nuigalway.ie Below we discuss their results. 
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Table 1. Log of observations of BF Eri 



Set 


Date 


HJD Start 


Instrument 


AA a 


A range 


Exp. Time 


Number 


Duration 






2450000+ 




(A) 


(A) 


(sec) 


of exps. 


(hours) 


2005-N1 


2005-Oct-30 


3673.982 


B&Ch 6 


6.0 


3680-6770 


300 


16 


1.44 


2005-N2 


2005-Oct-31 


3674.821 


B&Ch 


6.0 


6025-8000 


300 


40 


3.50 


2005-N3 


2005-Nov-Ol 


3675.821 


B&Ch 


2.1 


6145-7225 


300 


58 


5.25 


2006-Nov 


















2006-N1 


2006-Nov-22 


4061.884 


Echelle 


0.234 


3915-7105 


1200 


9 


3.13 


2006-N2 


2006-Nov-23 


4062.728 


Echelle 


0.234 


3915-7105 


1200 


17 


6.81 


2006-N3 


2006-Nov-24 


4063.710 


Echelle 


0.234 


3915-7105 


1200 


19 


7.08 


2006-N4 


2006-Nov-25 


4064.684 


Echelle 


0.234 


3915-7105 


1200 


3 


0.70 






4064.879 


Echelle 


0.234 


3915-7105 


1200 


3 


0.70 


2006-N5 


2006-Nov-26 


4065.698 


Echelle 


0.234 


3915-7105 


1200 


20 


7.26 


2006-N6 


2006-Nov-27 


4066.689 


Echelle 


0.234 


3915-7105 


1200 


14 


4.78 


2006-Dec 




















2006-Dec- 13 


4082.840 


Echelle 


0.234 


3915-7105 


1200 


2 


0.35 




2006-Dec- 14 


4083.850 


Echelle 


0.234 


3915-7105 


1200 


2 


0.35 




2006-Dec- 15 


4084.808 


Echelle 


0.234 


3915-7105 


1200 


2 


0.36 



a - AA is the FWHM spectral resolution 
6 - B&Ch - Boiler & Chivens spectrograph 



our observations and data reduction. The data analysis and the re- 
sults are presented in Section [3] Here we find the orbital period, 
measure the radial velocities of the WD and the secondary star, and 
determine the spectral type of the secondary and its rotational ve- 
locity. In Section [4] we discuss the results of Doppler tomography 
of the emission lines from BF Eri. The system parameters are ob- 
tained in Section[5] while a discussion and a summary are given in 
Sections [6] and|7] respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

The spectra presented here were obtained at the Observatorio As- 
tronomico Nacional (OAN SPM) in Mexico on the 2.1m telescope. 
The first observations, in 2005, were conducted with the Boiler & 
Chivens spectrograph, equipped with a 24 /im (1024 x 1024) SITe 
CCD. Observations were made during three consecutive nights of 
October 30-31 - November 1 in the wavelength range of 3680- 
6770 A, 6025-8000 A and 6145-7225 A respectively. Spectra were 
obtained in the first order of 400 line/mm (for the first two nights) 
and 1200 line/mm (for the third night) gratings with corresponding 
spectral resolutions of around 6.0 A and 2. 1 A respectively. A to- 
tal of 114 spectra were obtained with 300 s individual exposures, 
half of them on the third night. Cu-He-Ne-Ar lamp exposures were 
taken before and after the observations of the target for wavelength 
calibrations, and the standard spectr ophotometric stars FeigellO, 
HR3454 and G191-B2B Jokej|l990) were observed for flux cali- 
brations. 

In order to check the current photometric state of the object, 
during the first night we obtained a few photometric observations 
on an accompanying 1.5 m telescope at the same site. The V mag- 
nitude was about 15.0-15.2, indicating that BF Eri was definitely 
in the quiescent state. These observations showed that BF Eri is a 
long-period system and that the obtained data were not sufficient to 
determine the orbital period with good precision. 

Further observations were obtained during 6 consecutive 
nights of November 22-27, 2006 and 3 nights of December 13-15 , 
2006 with the Echelle spectrograph jLevine &Chakrabartvll 19951) 



attached to the same 2.1 m telescope. This instrument gives a res- 
olution of 0.234 A pixel - at Ha using the UCL camera and a 
CCD-Tek chip of 1024x1024 pixels with a 24 ^m 2 pixels size. 
The spectra cover 25 orders and span the spectral range 3915- 
7105 A. A total of 85 spectra were obtained in November and 6 
spectra in December with 1200 s individual exposures. In Novem- 
ber we also took spectra of the spectral type templates Gliese 4. 1 
(G5V), WO 9808 (G6V), NN4366 (G7V), Gliese 793.1 (G9V), 
Gliese 75 (K0V), Gliese 68 (K1V), Gliese 33 (K2V), Gliese 183 
(K3V), Gliese 53.1 (K4V), Gliese 69 (K5V), Gliese 169 (K7V), 
V Psc (Ml V) and Gliese 844 (M2V). All the nights of observations 
were photometric and the seeing ranged from 1 to 2 arcsec. TableQ] 
provides a journal of the BF Eri observations. 

The reduction procedure was performed using IRAF. Compar- 
ison spectra of Th-Ar lamps were used for the wavelength cali- 
bration. The absolute flux calibration of the spectra was achieved 
by taking nightly echellograms of the standard stars HD93521 and 
HR153. Though we used a wide slit (2") with seeing usually no- 
ticeably less than the slit width, this does not warrant excellent 
flux calibration, since only an average curve for atmospheric ex- 
tinction and a permanent E-W orientation of the slit were used. 
At the same time, due to an unexpectedly observed outburst of 
BF Eri in the November set we found it useful to obtain some 
photometric and colour information from our spectra. For this we 
followed an approach whic h we used for investigating an outburs t 
of the dwarf nova BZ UMa dNeustroev. Zharikov. & Michelll2006t) . 
Namely, we defined an internal photometric system comprising 
four colour bands u b' v r' centered at 4000 A, 4550 A, 5590 A, 
6380 A with widths of 50 A, 100 A, 120 A, 150 A respectively. The 
colour indices were calculated as C = —2.5 log(/i//2), where f\ 
and f% are the fluxes averaged across the corresponding bands. To 
check the stability and the flux calibration accuracy we have de- 
termined the u — b', b' — v' and v' — r' colours of the control 
star HR153 using our nightly spectra, and compared these with its 
published spectral energy distribution. The colours did not differ 
by more than 0.03 mag between any of our observations, and our 
colours were within 0.06 mag of the published spectra. 



Spectroscopy of the cataclysmic variable BF Eri 



3 



E 

jj 

'en 
E> 

CD 

o 

X 

3 




4000 4500 5000 5500 6000 6500 

Wavelength (A) 



7000 



7500 



8000 



Figure 1. The average spectra of BF Eri from the 2005 set, uncorrected for orbital motion. The lower spectrum is combined from the 2005-N1 and 2005-N2 
sets while the inset spectrum is an average of all spectra from the 2005-N3 set. 



To improve the confidence of the results presented in this pa- 
per we also acquired the night- and sef-averaged spectra obtained 
by means of co-adding of all spectra collected during each night 
and each set of observations respectively. Additionally we obtained 
the phase-averaged spectra for the 2005-N3 and 2006-Nov datasets. 
For this we have phased the individual spectra with the orbital pe- 
riod, derived in Section |33l and then co-added the spectra into 15 
and 17 separate phase bins respectively. In fact, only 11 bins were 
filled for the 2005-N3 set because these data were taken without 
complete orbital coverage. 



3 DATA ANALYSIS AND RESULTS 
3.1 Average spectra 

The mean spectra of BF Eri are shown in Fig.Q] These spectra are 
an average and a combination of all spectra from the 2005 set, un- 
corrected for orbital motion. During these observations the system 
appears to have been in a quiescent state (see Section [3~2t . The 
spectrum is typical for a dwarf nova. It is dominated by strong and 
broad emission lines of hydrogen and neutral helium. In addition 
to them He II A4686 is observed also but the C III/N III blend is not 
seen. All emission lines are single-peaked. The Balmer decrement 
is flat, indicating that the emission is optically thick as is normal 
in dwarf novae. The secondary star is clearly visible as absorption 
lines of the neutral metals Mg I, Fe I and Ca I, even in the average 
spectra, uncorrected for orbital motion. 

Much more detail can be seen in the Echelle spectra. Fig. [2] 
shows the selected averaged and normalised spectral orders from 
the 2006-Nov set of observations. We have selected to show those 
spectral orders which exhibit the strongest absorption-line features 
used for the following analysis. These plots were Doppler-corrected 
into the frame of the secondary star, employing the ephemerides 



from Table|5] The absorption spectrum is rich in features. In partic- 
ular we note the Mgl triplet at AA 5167-5173-5184 A (order 14), 
the Fel A5328 line (order 15), the Cal A6162 line (order 21), the 
Fel/Cal/Ball A6495 blend (order 23). Spectral order 4, although 
very noisy, indicates the presence of the Ca I A4226 line. We also 
additionally show order 23 which consists of the region around 
the strongest emission line Ha. While the signal-to-noise ratios of 
the absoption lines were increased after the Doppler correction, the 
emission line has been smeared out. To correct this, we also show, 
by the dashed line, the profile of Ha corrected for orbital motion of 
the WD. However note that even in this higher spectral resolution 
the emission line profile still appears to be single-peaked. 

Table |2]outlines fluxes, equivalent widths, velocity widths and 
relative intensities of the major emission lines measured from the 
averaged spectra of all the observational sets. 

3.2 Long term photometric and spectral changes 

At the beginning of the 2006-Nov set of observations we found 
the system to be appreciably brighter and bluer than in 2005, be- 
ing probably in an outburst. We have compared the fluxes in all 
four spectral passbands and found that the brightness of BF Eri in- 
creased by ~1.1 mag in it', ~0.8 mag in b', ~0.6 mag in v' and 
~0.4 mag in r' colours, relative to the 2005-N1 set (Fig. [3}. During 
the following nights we observed reddening of the flux distribu- 
tion and decreasing of the system brightness, which almost ceased 
towards the last two nights. Nonetheless BF Eri remained to be 
brighter and bluer than in 2005. However when we observed the 
system again two weeks later (the 2006-Dec set), we found it no- 
ticeably weaker than one year previously (Fig. [3}. 

The above-described photometric changes we observed 
throughout the outburst were accompanied by spectral changes 
which however were mostly quantitative rather than qualitative. 
There was no emergence of broad absorption troughs around the 
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Figure 2. The selected averaged and normalized spectral orders from the 2006-Nov set of Echelle observations of BF Eri, compared with several spectral 
standards. The spectra of BF Eri have been corrected for orbital motion of the secondary star. We also additionally show spectral order 23 which consists of the 
region around the strongest emission line Ha. While the signal-to-noise ratios of the absoption lines were increased after the Doppler correction, the emission 
line has been smeared out. To correct this we also show, by the dashed line, the profile of Ha corrected for orbital motion of the WD. The spectra of standards 
have been broadened by 78 kms" 1 to account for the rotational broadening of the lines in BF Eri and by 59 kms -1 to account for orbital smearing of the 
BF Eri spectra during exposure (see Section |3~71 and footnote |2]f or details). 



emission lines nor dr amatic changes in their profiles, as w e ob- 
served in BZ UMa dNeustroev. Zharikov. & Michell l2006h . In a 
quantitative sense the changes became apparent by the decreasing 
of FWHM, EW and flux of the emission lines during the first nights 
of the 2006-Nov set (Fig.[4] Table|2j- Also of particular interest is 
tracing the changes of the high excitation lines (such as the He II 
and C III/N III blend emissions) which are good tracers of irradia- 



tion. Unfortunately, due to the weakness of these lines and the poor 
sensitivity of the CCD at short wavelengths we were able only to 
detect the (probable) presence of He II A4686 in the outburst and 
no signs of the C III/N III Bowen blend. Thus, unlike many other 
dwarf novae which show strengthening of the He II and C III/N III 
line emissions during an outburst, these lines have remained very 
weak or even absent in the outburst spectra of BF Eri. 
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Table 2. Fluxes, Equivalent Widths (EW), Full Widths at Half Maximum (FWHM), Full Widths at Zero Intensity (FWZI) and Relative Intensities of the major 
emission lines 



Set 


Spectral 


Flux 


EW 


Relative 


FWHM 


FWZI 




line 


(xlO~ 14 ergs -1 cm -2 ) 


(A) 


Intensity 


(km s _1 ) 


(km s _1 ) 


2005 -JN 1 


Ha 


10.5 


32.3 


2.18 


1000 


4200 




KB 


16.0 


41.6 


2.41 


1320 


6600 






12.6 


39.6 


2.40 


1500 


4000 




He I A5875 


2.8 


7.5 


1.33 


1100 


4000 




He I A6678 


1.6 


5.1 


1.16 


1430 


3200 




He II A4686 


2.8 


7.7 


1.14 






2005-N2 


Ho 


11.3 


33.3 


2.46 


915 


3800 




He I A6678 


1.7 


5.3 


1.18 


1300 


3200 




He I A7065 


1.2 


3.9 


1.12 


1270 


3100 


2005-N3 


Ho 


12.7 


32.5 


2.50 


915 


3850 




He I A6678 


1.8 


4.9 


1.18 


1250 


2800 




He I A7065 


1.3 


4.0 


1.13 


1170 


3100 


2006-N1 


Ha 


17.6 


22.5 


2.30 


710 


2300 




H/3 


17.4 


12.2 


1.89 


800 


2350 




He I A5875 


3.6 


3.6 


1.22 


870 




2006-N2 


Ha 


9.8 


15.7 


1.97 


670 


1900 




H/3 


7.4 


6.5 


1.54 


710 


1650 




He I A5875 


1.7 


2.0 


1.17 


670 




2006-N3 


Ha 


11.0 


18.2 


1.99 


670 


2900 




H/3 


7.1 


5.7 


1.45 


650 


2300 




He I A5875 


2.0 


2.2 


1.18 


740 




2006-N4 


Ha 


8.9 


15.5 


1.92 


660 


2350 




II . 


7.1 


6.9 


1.49 


800 


1800 




He I A5875 


2.0 


2.6 


1.16 


690 




2006-N5 


Ha 


9.3 


16.3 


1.92 


650 


2500 




H/3 


9.9 


10.7 


1.66 


1050 


2500 




He I A5875 


2.5 


3.5 


1.22 


890 




2006-N6 


Ha 


12.3 


23.0 


2.12 


760 


3000 




H/3 


8.2 


8.5 


1.55 


960 


2350 




He I A5875 


2.3 


3.2 


1.18 


820 




2006-Nov 


Ha 


10.8 


17.0 


2.02 


670 


2250 




H/3 


8.3 


7.5 


1.51 


830 


2450 




H 7 


7.3 


5.3 


1.42 


920 






He I A5875 


2.2 


2.6 


1.17 


790 


1600 




He II A4686 


1.6 


1.2 


1.10 






2006-Dec 


Ha 


17.4 


26.0 


2.46 


730 


2550 




H/3 


19.3 


25.0 


2.35 


950 


2700 




He I A5875 


5.6 


8.2 


1.35 


1430 


2000 



3.3 Orbital period determination 

In order to determine the orbital period of BF Eri we measured 
the radial velocity variations of Ha by convolving the observed 
line profiles with a single Gaussian of FWHM = 600 km s . 
The spectra were continuum- normalized prior to this analysis. 
The data were studied for period icities using the standard Lomb- 
Scargle power spectrum analysis dLombll 19761 ; [Scargle 1982). Due 
to the longer observing span and more homogeneous nature of 
the data we initially studied the observations from the 2006 set. 
The resulting periodogram shows a strong signal at the 3.6940 
day -1 frequency (Fig. [5), which corresponds to the orbital period 
of P rt = 0.2707 ± 0.0008 days. 



A more accurate value of the orbital period was deter- 
mined Jrojn_Jhe_rambmed_date ; _J^e applied the CLEAN proce- 
dure ( Ro berts. Lehar. &Drehej[l987h to sort out the alias peri- 
ods resulting from the uneven data sampling and daily gaps and 
obtained a strong peak in the power spectrum at the 3.69165 
day -1 frequency, corresponding to the orbital period of P or t = 
0.270 881 ±0.000003 day s (Fig.g), consistent with the determina- 
tion of lSheets etall ( 120071) . 

3.4 The radial velocity of the white dwarf 

In CVs the most reliable parts of the emission line profile for de- 
riving the radial velocity curve are the extreme wings. They are 
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Figure 3. The u' — b', b' — v' and v' — r' colour indices observed in 
BF Eri. Symbols represent data obtained from the m'gfa-averaged spectra. 
The v' light curve is also shown (dots, right Y). 
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Figure 4. Equivalent Width (EW), line fluxes and Full Widths at Half Max- 
imum (FWHM) variations of the Ha, H/3 and He I A5875 emission lines 
observed in BF Eri. 



presumably formed in the inner parts of the accretion disc and 
therefore should represent the motion of the WD with the high- 
est reliability. We measured t he radial velocities using th e double- 
Gaussian method described bv lSchneider & Young| ( 1980) and later 
refined bv lShaftedfl983h . This method consists of convolving each 
spectrum with a pair of Gaussians of width a whose centres have 
a separation of A. The position at which the intensities through 
the two Gaussians become equal is a measure of the wavelength 
of the emission line. The measured velocities will depend on the 
choice of a and A, and by varying A different parts of the lines 
can be sampled. The width of the Gaussians a is typically set by 
the resolution of the data. In order to test for consistency in the de- 
rived velocities and the zero phase, we separately used the emission 
lines Ha and H/3 in the 2006-Nov phase-averaged spectra, and the 
emission line Ha in the phase-averaged spectra from the 2005-N3 




Frequency (day" ) 

Figure 5. Lomb-Scargle power spectrum of the Ha radial velocity varia- 
tions from the 2006 dataset. The inset shown is blow-up of the most promi- 
nent peak of the CLEANed power spectrum of the combined data. 



set. All measurements were made using a Gaussian FWHM of 100 
km s -1 (and also of 50 km s _1 for the 2006 dataset) and differ- 
ent values of the Gaussian separation A ranging from 300 km s _1 
to 2000 kms -1 in step s of 50 kms" 1 , following the techni que of 
"diagnostic diagrams" dShafter. Szkodv, &"T horstensen 198JJ). 

For each value of A we made a non-linear least-squares fit of 
the derived velocities to sinusoids of the form 



V(4>, A) = 7 (A) - iv"i(A) sin [2tt (<j> - <t> (A))] 



(1) 



where 7 is the systemic velocity, K\ is the semi-amplitude, <j>o is 
the phase of inferior conjunction of the secondary star and <f> is 
the phase calculated relative to epoch To = 2453675.2670 for the 
2005-N3 set and T = 2454061.8112 for the 2006-Nov set. 

The resulting "diagnostic diagram" for Ha in the 2006 dataset, 
with a — 100 km s _ , is shown in Fig. [6] The diagram shows 
the variations of K\, a(K\)/Ki (the fractional error in Ki), 7 
and 4> with A dShafter. Szkodv. & Thorstensen|[l986h . The dia- 
grams for other lines and a look the same. To derive the orbital el- 
ements of the line wings we took the values that correspond to the 
largest separation just b efore a(K\)fK\ shows a sharp increase 
dShafter & Szkodvll 19841) . We find the maximum useful separation 
to be A ~ 1100-1250 kms -1 . In Fig. [6]it appears that A can 
be increased to ~ 1 150—1250 km s _1 before a(K)/K begins to in- 
crease. We also note that all the calculated orbital elements are very 
stable, even for the smaller separation A, hence there are no diffi- 
culties in choosing their values. This statement is also true for other 
"diagnostic diagrams". We have obtained very consistent results for 
both the radial velocity semi-amplitudes and the 7-velocities for all 
of them, and we adopt the mean values K\ — 74 ± 3 km s -1 and 
7 = — 91 ± 7 km s _1 . The measured parameters of the best fitting 
radial velocity curves are summarized in Tables [3] and [5] In Fig. [7] 
we show the radial velocity curves for the Ha emission line. 

The derived value of the radial velocity semi-amplitude K \ 
is highly inconsistent with the one found by ISheets et ail d2007l) . 
We are unable to explain why their result is so different, as they 
did not present their analysis in detail. On the other hand, we have 
no reason to doubt our values. A possible reason might be poor 
spectral resolution for their data and as a consequence the wider 
Gaussians used in the double-Gaussian metho d. According to the 
"diagnostic diagram" (Fig. [6j we presume that lSheets et aljd2007t) 
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Figure 6. The diagnostic diagram for the Ha 2006-Nov data, showing the 
response of the fitted orbital elements to the choice of the double-gaussian 
separation. The best fit is reached with the gaussian separation of 1 150— 
1250 km s -1 . 



Table 3. Elements of the radial velocity curves of BF Eri derived from emis- 
sion lines 



Emission line & 
dataset 


Ki 
(kmr 1 ) 


7-velocity a 
(km s -1 ) 


0o 6 


Ha - 2005 


74±4 


-98±4 


0.013±0.011 


Ha - 2006 


74±4 


-66±5 


0.012±0.015 


H/3 - 2006 


77±10 


-73±12 


0.013±0.016 


Mean 


74±3 


-91±7 


0.0127±0.0003 



a The measured 7-velocities are heliocentric. The mean value was obtained 
after correction for the solar motion. 

b Phases were calculated relative to epoch To = 2453675.2670 for the 
2005 set and T = 2454061.8112 for the 2006 set. 

examined the lower- velocity parts of the Ha line and consequently 
obtained a higher value of K\ . 



3.5 The radial velocity of the secondary star 

The absorption lines from the secondary star in BF Eri are visi- 
ble quite clearly, even in the individual spectra of the 2005-N3 set, 
and more easily in the phase-folded spectra of both the 2005-N3 
and 2006-Nov sets. In order to obtain radial velocities, we used the 
following interactive method. 



(i) The first step was to measure the radial velocity variations of 
the absorption line Cal A6162 (and the Fel/Cal/Bal A6495 blend 
for comparison purpose) in the phase-folded spectra of the 2005- 
N3 set. The velocities were measured by convolving the observed 
line profiles with a single Gaussian. The resulting radial velocity 
curves were then fitted with a sinusoid of the form 



V(<j>) = 7 - K 2 sin [27r ( 



(2) 



and preliminary values of the radial velocity semi-amplitude K%, 
the systemic velocity 7 and the phase zero-point 0o were calcu- 
lated. The obtained parameters are K2 = 183 ± 6, 7 = —94 ± 5 
and <f> = 0.561 ± 0.007. 

(ii) Each spectrum in both the 2005-N3 and 2006-Nov sets was 
then shifted to correct for the orbital motion of the donor star, and 
the results averaged. These averaged spectra of BF Eri in the rest 
frame of the secondary star allowed us to estimate a preliminary 
spectral type of the donor star (see Section |3~6l . The averaged spec- 
tra of BF Eri were then cross-correlated with the velocity standards 
of the obtained spectral type and the relative shifts computed and 
applied. 

Thus, we obtained the cross-correlation template spectra which 
were used for measuring the absorption-line velocities. This was 
preferable to cross-correlating with the standard stars for the ob- 
vious reason that it maximizes the similarity between the template 
and the individual spectra to be cross-correlated, and avoids intro- 
ducing errors caused by spectral-type mismatch between the stan- 
dard star and the secondary star in BF Eri, as well as by the in- 
correctly determined rotational broadening of the absorption lines. 

(iii) Next, the phase-folded BF Eri spectra were cross-correlated 
with the templates. For the 2005-N3 set, the wavelength region 
6150-6515 A was used. To avoid the influence of the night-sky 
lines A6300 A and A6363 A, some portions of the spectral range 
around these lines were blanked. For the Echelle 2006-Nov set, the 
procedure was carried out using orders 4, 14, 15, 16 and 21 in the 
wavelength regions 4220-4235 A, 5160-5235 A, 5260-5400 A, 
5360-5500 A and 6100-6270 A respectively (Fig.|2j. The solutions 
obtained by fitting the measured radial velocities with sinusoid $2^ 
are very similar for both sets and all orders. 

(iv) Steps |(ii)1 and |(iii)1 were then repeated and the final results 
(with almost no changes from the previous step) were obtained, 
so we adopt the averaged values of K2 = 182.5 ± 0.9 kms -1 
and 7 = —93.6 ± 0.4 kms -1 . These values are consistent with 
ISheets et all d2007t) . We also note that the difference between the 
phase zero-points obtained from the emission and absorption lines 
is very close to 0.5, as it must be if the derived velocities from 
those lines trace the components' motion. In Fig. UJwe show the 
measured radial velocities together with their sinusoidal fit. The 
measured parameters of the best fitting radial velocity curves are 
summarized in Tables|4]and|5] 

(v) Finally, the step |(iii)| was repeated again, using the stan- 
dard templates of the best-fitting spectral type K3V broadened to 
the predicted rotational velocity of the mass donor V ro t sin i=78 
kms -1 (see Section |3~71 and footnote [2] for details). No significant 
difference in the measured parameters was found. 



3.6 Spectral type of the secondary star 

The best absorption lines for spectral classification in normal 
stars, which are independent of chemical abundances are the line 
ratios of the Fel lines A4250, A4260, A4271 to the CrI lines 
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Figure 7. Radial velocities of the white dwarf and the secondary star in BF Eri folded on the ephemeris from Table[5] Observations were performed in 2005 
November (left panel) and 2006 November (right panel). Two cycles are shown for clarity. 

Table 4. Elements of the radial velocity curves of BF Eri derived from absorption lines by means of cross-correlation. 



Dataset 


Spectral 


Wavelength 


K 2 


7-velocity a 






Order 


region (A) 


(kms" 1 ) 


(kms" 1 ) 




2005 




6150-6515 


185±5 


-94±4 


0.500±0.005 


2006 


Order 04 


4220-4235 


180±4 


-80±3 


0.502±0.003 


2006 


Order 14 


5160-5235 


187±3 


-83±2 


0.499±0.002 


2006 


Order 15 


5260-5400 


183±2 


-81±1 


0.499±0.001 


2006 


Order 16 


5360-5500 


182±2 


-81±2 


0.500±0.002 


2006 


Order 21 


6100-6270 


181±2 


-82±1 


0.499±0.001 


Mean 






182.5±0.9 


-93.6±0.4 


0.500±0.001 



a The measured 7-velocities are heliocentric. The mean value was obtained after correction for the solar motion. 
b Phases were calculated relative to epoch T = 2453675.2670 for the 2005 set and T = 2454061.8112 for the 2006 set. 



A4254 and A4274 as well as the strength of th e Ca I A4226 line 
dKeenan & McNeill 19761 : lEchevarria et al.l2007l) . The Cr I lines in- 
crease in strength with respect to the Fe I lines with spectral type. 
On the other hand the intensity of the Cal line also increases 
steadily as the spectral type advances (in fact this line has a satura- 
tion effect for late K and M stars). Unfortunately, due to the poor 
sensitivity of the CCD at short wavelengths the data were too noisy 
in this wavelength region, so we were able to use only the width 
of the Ca I A4226 line as a spectral indicator. In accordance with 
this a spectral type of the secondary star in BF Eri can be initially 
estimated as K3-K4. 

To verify this result, we also analyzed the appearance of the 
absorption lines in those spectral orders where they are visible bet- 
ter. Purely by inspection, we have found several reliable indicators 
to be useful in constraining the spectral type. Among them are pairs 
of the lines/blends Fel A5227.2 and Crl/Fel AA5206+5208 (order 
14), Fel A5434.5 and Fel A5424.1 (order 16), Cal A6122.2 and 
Fel A6137 (order 21), Cal A6439.1 and Fe I AA6400.0+6400.3 (or- 
der 22). The relative depths of all these couples change rapidly over 
the spectral type range we are interested in here. Analysis of these 
absorption features supports a K2-K3 spectral classification for the 
donor star in BF Eri. 

Taking into account all the estimations we adopt the final value 
for the spectral type of the secondary to be K3±0.5, consistent with 
the determination of Shee ts etal] d2007l) . 



3.7 Rotational velocity of the secondary star 

The absorption lines visible in the averaged and orbitally corrected 
spectra are obviously broadened due to orbital smearing during ex- 
posure and rotation of the secondary star. The observed rotational 
velocity of the secondary could provide an independent determina- 
tion of the mass ratio of the binary. We have attempted to estimate 
this parameter by artificially broadening the template star spectra 
which are assumed to have low V ro t sin i, and fitting them to the 
orbital velocity corrected BF Eri spectra until the lowest residual 
was obtained. We have restricted this analysis to the spectra ac- 
quired with the Echelle spectrograph because of their higher spec- 
tral resolution. We have selected from them those spectral orders 
which exhibit the greater set of the strong absorption-line features 
(orders 14, 15, 16 and 21). Al l details of the used V ro t sin i m easur- 
ing technique can be found in Smith. Dhillon. & Marsh 1 19981) [see 



also lMarsh. Robinson, & Woo j Jl994h and lNorth et all feOOCh l. 

Unfortunately, despite the efforts expended we have failed to 
derive consistent value for the rotational velocity. Instead we ob- 
tain a quite broad range of V ro t sin i values of 70-100 km s~ , that 
is unacceptable for the determination of the system parameter^ 

2 See corresponding discussion in Section |6~T1 Note that using the system 
parameters of BF Eri, obtained in Section|5](Table[5j, we get the estimated 
value of Vrot sin i to be ~ 78 km s , consistent with observations. 
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We have additionally examined the spectra, varying the wavelength 
ranges and even cutting the separate absorption lines, but were un- 
able to get a more certain value of V ro t sin i. We explain our fail- 
ure by the insufficient quality of the used spectra as even the aver- 
aged spectra are quite noisy. We still need more observations with 
high spectral resolution and better time resolution to fully solve this 
problem. 

Although we could not determine the rotational velocity, we 
were able to estimate the fractional contribution of the secondary 
star to the total light. We have found that the secondary contributes 
about 1 1 per cent of the total light in order 14. 



4 DOPPLER TOMOGRAPHY 

The orbital variation of the emission line profiles indicates a non- 
uniform structure for the accretion disc. In order to study the emis- 
sion structure of BF Eri we have used Do ppler tomography. Ful l 
tec hnical details o f the method are given by Mar sh & Home! Jl988h 
and|Marsh (2001). Examples of the applicat ion of Doppler tomog- 
raphy to real data are given bv lMarsfi ( feOOll) . 

Figure [8] shows the tomograms of the Ha and He I A6678 
emissions from the 2005-N3 set and of the Ha, H/3, H7 and 
He I A5876 emissions from the 2006- Nov se t of ob servations, com- 
puted using the code developed by ISpruitl Il998h . These figures 
also show trailed spectra in phase space and their corresponding 
reconstructed counterparts. A help in interpreting Doppler maps 
are additional inserted plots which mark the positions of the WD 
(lower cross), the center of mass of the binary (middle cross) and 
the Roche lobe of the secondary star (upper bubble with the cross). 
The Roche lobe of the secondary has been plotted using the system 
parameters, derived in Section[5] 

Due to the non-double-peaked emission line profiles of 
BF Eri, we did not expect a Doppler map to have an annulus of 
emission centered on the velocity of the WD, and the observed 
tomograms do not show it. Instead of this all the maps display a 
similar, very nonuniform distribution of the emission, the bulk of 
which is located in the lower-left quadrant of the tomograms, at 
lower velocities than the predicted outer disc velocities for an ac- 
cretion disc with a radius Rl 1 ■ The appearance of the 2005-N3 Ha 
map is similar to the hydrogen line maps from the 2006-Nov set 
but a most prominent contribution of the emission here is from the 
area around the center of mass of the binary system. Note also that 
BF Eri does not show any emission on the hemisphere of the donor 
star facing the WD and boundary layer. 

The interpretation of the emission structure in BF Eri is am- 
biguous. All the detected emission sources are located far from 
the region of interaction between the stream and the disc parti- 
cles. This 'reversed bright spot' phenomenon can perhaps be ex- 
plaine d by a gas stream which passes above the disc and hits its 
back dHeilier & Robinsonll 19941) , or alternatively, by the disc thick- 
ening in resonating locations. 



5 THE BINARY SYSTEM PARAMETERS 

It is impossible to determine accurate system parameters for a non- 
eclipsing binary system like BF Eri. However, we can constrain 
them, using the measurements of the radial velocities of the primary 
and secondary stars K\ and K2 in conjunction with the derived 
orbital period P or b and the spectral class of the secondary. 



Combining the values of K\, K% and P we find the masses of 
each component of the system 

Mr sin 3 i = + K # = 0.34 ± 0.01M a , (3) 

2ttG 



M 2 sin i = 



PK 1 {K 1 +K 2 f 
2-kG 



0.14±0.01M, 



■ 



and the projected binary separation 



a sin % — 



P{K 1 +K 2 ) 



2tt 



= 1.37±0.02i? G 



(4) 



(5) 



The inclination angle i of the system is an unknown param- 
eter but for a CV it can be restricted by using some reasonable 
assumptions. The published photometry of BF Eri is extensive 
enough to rule out any significant eclipse, which constrains the 
inclination i to be less than about 70° to avoid an obvious par- 
tial eclipse of the disc. At i = 70° the minimum WD mass is 
around 0.41 Mq. The Chandrasekhar mass limit determines the 
lowest limit for the inclination angle to be about 38°. Thus, we 
have now the pure dynamic solution for the parameters of BF Eri: 
q = Ms /Mi = K1/K2 = 0.41 ± 0.02, Mi=0.41-1.44 M©, 
M 2 =0.17-0.59 Mq, i=38°-70°, a=l. 46-2.23 R e . 

At this point we would like to take note that from general theo- 
retical considerations the donor star in a long-period CV should be 
more massive than in a short-period one. In practice, using any of 
the rec ently obtained empirical mass-period relations (see, for ex- 
ample. lWarneJll995l : ISmith & Dhillonlll998l ; |Patterson et alj|2005l) 
we obtain the mass for the secondary 0.67-0.72 Mq (Fig. [9] bot- 
tom panel). This is even more than our upper bound for the donor 
mass, leading us to the conclusion that the WD in BF Eri is very 
close to the Chandrasekhar limit. 

Also we can estimate the secondary's mass, reasonably assum- 
ing that the secondary star fills its Roche lobe. The relative size of 
the donor star is constrained by Roche geometr y and can be esti- 
mated using Eggleton's formula (Eggleton 1983) 



R.2 



0.49g 



2/3 



(6) 



a 0.6(7 2 / 3 + ln(l + q 1 ' 3 ) ' 

which gives the volume-equivalent radius of the Roche lobe to bet- 
ter than 1 per cent. 

In Fig. [9] (upper panel) we show the Roche lobe size of the 
secondary as a function of the binary separation, and also differ- 
ent theoretical and empirical Mass-Radius relations. One can see 
that in order to fill its Roche lobe, the secondary must be relatively 
massive (0.7-0.8 Mq) if it is on the main sequence. But the mass of 
the primary in this case must exceed the Chandrasekhar limit, con- 
tradicting observations. A very massive white dwarf might just be 
accommodated again. However the secondary of the corresponding 
mass of 0.59 Mq must be significantly larger than a main-sequence 
star of the same mass. It is now well known that all CVs have sec- 
ondaries slightly larger (~ 1 5 per cent) than the zero-age main- 
sequence (ZAMS) single-star jPatterson et alJl2005l;lKnigge1l2006l) 
but the secondary in BF Eri is even greater by another 15 per cent. 

Additional constrains on the mass of the secondary can be ob- 
tained using another independently obtained parameter - its spec- 
tral type. Actually, this method should be used with great caution as 
it is understood that there is a big range in mass for a given spectral 
type JSmith & Dh illon 1998). Ho wever, using the latest g eneration 
of low-mass star models, iKolb. King. & Baraffd ( l200ll) explored 
theoretically to what extent the spectral type provides a reasonable 
estimate of the donor mass. They concluded that the spectral type 
should be a good indicator of the donor mass, and gave lower and 
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Figure 8. Doppler tomography for the Ho and He I A6678 emission lines from the 2005-N3 set and for the He I A5876 line from the 2006-Nov set (in the upper 
half of Figure), and for Ha, H/3 and H7 from the 2006-Nov set of observations (in the bottom half of Figure). For each line the observed and reconstructed 
trailed spectra (bottom) and corresponding Doppler maps (top) are shown. 
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Table 5. Orbital and system parameters for BF Eridani. 



Parameter 


Measured 


Monte Carlo 




Value 


Value 


Porb (d) 


0.270881±3 x 10~ 6 




T (+2450000) 


4061.81 12±0.0003 




K\ (kmr 1 ) 


74±3 




K 2 (kms- 1 ) 


182.5±0.9 




7 (km s~ ) 


-93.6±0.4 




Spectral Type 






of secondary 


K3±0.5 




d (pc) 


700±200 




Mi/M 


1.23-1.44 


1.28±0.05 


M 2 /M Q 


0.50-0.59 


0.52±0.01 


q = Afa/Mi 


0.41±0.02 


0.41±0.01 


i 


38°^10.5° 


40° ±1 


a/R Q 


2.11-2.23 


2.14±0.02 



upper limits to M2 for a given spectral type (see Table 1 in their pa- 
per). For the secondary star in BF Eri we have the range of masses 
of ~ 0.5 — O.8M0 . The upper limit here is for the mass of a ZAMS 
star with the same spectral type as the donor while the lower limit is 
for evolved mai n sequence stars under some extreme assumptions 
as discussed bv lKolb. King, & Baraffel J 200 lh . Thus, this leads us 
to conclude that the secondary in BF Eri is an evolved star of the 
mass of ~ 0.50 — O.59M0. This in turn explains why the size 
of the secondary is so big for its mass. As Baraffe & Kolb (2000) 
showed, the size of the donor star of any given mass is larger for 
more evolved sequences. 

Thus, from the above analysis based on strong assumptions, 
we can restrict the system parameters of BF Eri in the following 
ranges: q=0.41, Mi=l, 23-1.44 Mq, M 2 =0.50-0.59 Mq, i=38°- 
40.5°, a=2. 11-2.23 Rq. Finally, in order to esimate the "best- 
fitting" values of the component masses and their errors, we use d 
a Monte Carlo approach similar to lHorne. Welsh. & W ade ( 1993). 
The Monte Carlo simulation takes 10 6 sample values of K\ and 
K2, treating each as being normally distributed about their mea- 
sured values with standard deviations equal to the errors on the 
measurements. The values of the inclination i are sampled at ran- 
dom from an uniform distribution over the interval [0, imax], where 
imax = 70°. Using the relations and (|4j, we then calculate the 
masses of the components. If Mi exceeds 1.44 Mq, or M2 is less 
than 0.50 Mq, the outcome is rejected. The ensuing distribution for 
accepted outcomes is used to calculate the component masses. The 
values of all the system parameters deduced from the Monte Carlo 
simulation are listed in Table[5] Each value corresponds to the peak 
and standard deviation of the relevant distribution. 



6 DISCUSSION 

6.1 System parameters 

We would like to make some final remarks on the system param- 
eters of BF Eri, derived in Section [5] These parameters were de- 
termined using a traditional spectroscopic solution based on the 
derived radial velocity semi-amplitudes. In order to calculate the 
masses of the binary we have assumed that the semi-amplitudes of 
the measured radial velocities represent the true orbital motion of 
the stars. In CVs, however, this may not be an accurate assump- 
tion, as the emission lines arising from the disc may suffer several 




Figure 9. Constraints on the system parameters of BF Eri. In the bottom 
panel, we show the Mi - M2 diagram. The diagonal solid line correspond 
to the derived mass ratio q=0.4 1 . The dots along this line correspond to dif- 
ferent mass solutions for specific values of the inclination angle, labelled 
nearby to those dots. The two vertical areas (shown in grey) mark the re- 
gions of impossible values, due to lack of eclipses (on the left) and the 
Chandrasekhar mass limit (on the right). The two horizontal lines s how the 
predic ted mass of the seco n dary f rom the M^-Porb relations by Warner 
119951) and IPattersonetaP <2005l) (the dashed and dotted lines, respec- 
tively). The green vertical line pattern circumscribes the range of masses of 
the evolved main sequence st ar of the spectral type K 3 under some extreme 
assumptions as discussed by iKolb. King. & Baraffe 12001). In the upper 
panel we show, by the thick solid line, the Roche lobe size of the secondary 
flasa function of the mass of the WD Mi, for the mass ratio g=0.41 and 
P or .(,=0. 270881 days. The lower lines correspond to different theoretical 
and empirical Mass- Radius relatio ns. The thick dashed line is the theoreti- 
cal M-R relation of lBaraffe et~al Jl998|) The th in dashed and dash - dotted 
lines are the empirical relations of Knigge ( 2006) and Patterson et al. ( 2005) 
respectively. The short-dotted line is the M-R relation based on photomet- 
ric, geometrical, and absolute elements for 112 eclipsing detached binary 
systems with both components on the main sequ ence and with kno wn pho- 
tometric and spectroscopic orbital elements IGorda & Svechnikov 1998). 



asymmetric distortions. Likewise, absorption line profiles may be 
subjected to irradiation or hot-spot contaminations in the surface of 
the secondary, distorting therefore the true value of K%, 

Taking these potential sources of errors into consideration we, 
however, believe we could avoid them. Indeed, though those prob- 
lems were severe in the early determinations of the radial velocities, 
nevertheless modern high resolution spectroscopy has greatly im- 
proved the methods in detecting asy mmetric disto rtions to provide 
more reliable radial velocity values dWarneJI 19951) . In our study of 
BF Eri, the value of K2 is beyond question. Scattering of its values 
derived with the use of the different data, is very small. Doppler 
tomography does not show any emission from the secondary star 
that can distort the absorption line profiles from the latter, neither 
in quiescence nor in the outburst. Furthermore, a non-uniform ab- 
sorption distribution across the surface of the secondary would re- 
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suit in a non-sinusoidal radial velocity curve that is not observed in 
BF Eri (Fig. [7J. All the obtained values of K\ are also consistent 
with each other. Though the detected asymmetric emission struc- 
ture of BF Eri (particularly, the spot in the lower-left of the Doppler 
maps) may potentially influence the velocity determination, we be- 
lieve we could avoid this as that strong emission spot is situated 
well inside of the chosen Gaussian separation. The correct phasing 
of the radial velocity curve also strengthens our confidence. 

Alternatively, one may also calculate the mass ratio of BF Eri 
using the observed rotational velocity of the secondary star, inde- 
pendently of K\ . Unfortunately, we were unable to derive a consis- 
tent value of the rotational velocity, getting instead a broad range 
of Vrot sin i values. This gives us the q values in the range of 0.33- 
0.62 that is consistent with the former solution, but also very un- 
certain, so we decided to reject it. 

Thus, we have shown that BF Eri contains a massive WD and 
an evolved secondary. This places this system among a small group 
of CVs with a high-mass WD (such as RU Peg and U Sco). Such 
systems are specially interesting sin ce they are considered as pos- 
sible Type la supernova progenitors jParthasarathv et al]|2007l and 
references therein). The identity of the progenitor has become im- 
portant for cosmology, since Type la supernovae are used as stan- 
dard candles and their peak luminosity depends on the nature of 
the progenitor. The exploding star is now thought to be a carbon- 
oxygen WD that accretes mass from a binary companion until it 
approaches the Chandrasekhar limit, ignites carbon under electron- 
degenerate conditions, undergoes a thermonuclear instability, and 
disrupts completely. Further studies of the systems with massive 
WDs similar to BF Eri are needed to understand if they evolve to 
SNe la. 



6.2 7-velocity, proper motion and distance 

A decade ago Ivan Paradiis, Augusteiin. & Stehlel dl996l) drew at- 
tention to the distribution of space velo cities, as a means of prob- 
ing the age profile of the CV population. iKolb & Stehlel d 19961) de- 
rived the theoretically expected distribution of 7-velocities and the 
dispersion of 7 as a function of orbital period. In particular, they 
showed that the CVs having periods longer than the upper limit of 
the period gap (P or b ^3 h) should be younger (^1.5 Gyr), and 
therefore have a smaller line-of-sight velocity dispersion according 
to the empirical age-velocity dispersion relation (predicted value 
^15 kms -1 ). Conversely, those CVs with orbital periods shorter 
than the lower limit of the period gap, should be older Qs 3-4 Gyr) 
and show a larger velocity dispersion (predicted as ~30 kms -1 ). 
Paying a s pecial attention to determine accurate absolute systemic 
velocities. iNorth et alj J2OO2T) found its values for four long-period 
DNs. They obtained an average 7-velocity of ~5 km s _1 with a dis- 
persion of ~8 km s _1 indicating that the above postulate is indeed 
true. 

In our study of BF Eri, we have mostly followed the method of 
INorth et al.1 d2002h and found the absolute 7-velocity, after correc- 
tion for the solar motion, to be —94 km s _1 , referred to the dynam- 
ical local standard of rest (LSR). This value is very far away from 
any expectations of the theory (as set out bv lKolb&Stehie1ll996h . 
In principle, such a kick velocity could be caused by an event in 
which a significant fraction of the total mass of the system is vio- 
lently, asymmetrically ejected. In the case of CVs, the progenitors 
indeed eject most of their mass during the common-envelope pe- 
riod, but it is supposed that this process doe s not alter the space 
velocity of the system. Kolb & Stehle dl996l) note, however, that 



this velocity may be affected by the cumulative effect of repeated 
nova eruptions with asymmetric envelope ejection. 

It is worthy to mention that BF Eri shows not only 
the high 7-velocity but also has a substantial proper motion: 
{ix x ,u,y) = (+34,-97 ) mas yr~ 1 dKlemola. Hanson. & Jonesl 
2000; Ha nson etalj|2004l) . In order to convert this value to the 
space velocity we need to know the distance to the system, which 
we can evaluate here by oblique methods only. 

One of these estimations is based on th e statistical period- 
luminosity-colours relation of lAk et all d2007l) . We have taken the 
JHK magnitudes of BF Eri from the 2MASS (Two Micron Al l 
Sky Survey) observations dHoard et al.l2002MSkrutskie et al.l2006h 
and found the distance to the system, from the dependence of the 
absolute magnitude Mj on the orbital period P or t and colours 
(J - H) and (H - K s ), to be 700 ± 200 pc. 

Alternatively, the distances can be estimated based on a semi- 
empirical donor sequence for CVs of iKn igge (2006). Unfortu- 
nately, Knigge's donor sequence is intended to describe the unique 
evolution track followed by unevolved secondaries, and he limited 
it to P or b < 6.2 h. BF Eri is a system with a longer period and 
likely an evolved secondary. However, we have slightly extrapo- 
lated Knigge's sequence and found from it and the JHK magni- 
tudes the lower limit of the distance d to be ~ 500 pc. Applying the 
empirical offsets to the donor absolute magnitudes (see, however, 
comments by Knigge on this) we have obtained d to be ~ 900 pc. 

Finaly, the distances can also be derived from the absolute 
magnit ude at outburs t My ( max) versus the orbital period rela- 
tion of lHarrison et al.l d2004l) . if one intends to consider the flares 
of BF Eri as DN outbursts (see corresponding discussion in Sec- 
tion |63). Analysis of the AAVSO database of BF Eri's observa- 
tions reveals 8 events brighter than 13 mag with an averaged value 
of 12.5 ± 0.2 mag. This yields the distance to be 650 ± 60 pc. 

Fortunately, all these independent estimations of the distance 
are consistent with each other, leading us to a final value of 
700 ± 200 pc. This corresponds to a transverse velocity vt = 
340 ± 100 kms" 1 or a space velocity of ~350 kms -1 with re- 
spe ct to the LSR , that i s again unusually high. Based on kinemat- 
ics, ISheets et al.l d2007t) surmised that BF Eri may belong to the 
halo population. Our Galactic components of the space velocity 
of the system (U, V, W) = (-2 20, -280, 4) kms'Ve slightly 
different from those obtained by ISheets et alj but also support a 
system motion being essentially in the Galactic plane and lagging 
far behind the rotation of the Galactic disc. In this connection the 
referee has pointed out that the high space velocity may be sim- 
ply a sign of great age and not a result of any asymmetric mass 
ejections, that would be consistent with the secondary star having 
started to evol ve prior to the time t he system came into contact. 
Indeed, though lKolb & Stehld Jl99^) do mention asymmetric nova 
explosions as a possible mechanism for increasing the dispersion in 
the 7-velocity distribution, but they give no numbers, and it is not 
yet clear whether it is really plausible that the effect could produce 
a 7-velocity as large as the one we found. 

However, if such a high space velocity of BF Eri was indeed 
caused by repeated nova eruptions then it might be worth searching 
for its past (recurrent) nova outbursts in the astronomical archive. 
Although the system is known only since 1940 and during the fol- 
lowing 50 years there were no observations, some collections of 
photographic plates cover over a century of data. It is straightfor- 
ward to look at plates taken over the years to see if past eruptions 
have been recorded. 
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Figure 10. Left: The AAVSO light curve of BF Eri. Blue stars represent the data labelled in the AAVSO database as V-magnitude, black filled dots - as 
Visual-magnitude, open dots - as inaccurate Visual-magnitude, green triangles - as the "Fainter Than" Visual-data. Right: Histogram of the reliable AAVSO 
data for BF Eri. 



6.3 The outburst behaviour and the nature of BF Eri 

The current classification of BF E ri as a DN is based on ob- 
servations by Watana be and Kato dWatanabej[l999l : lKatolll999l : 
iKato & Uemunfeood) who detected "outbursts" in the light curve 
of the system. We note that all those outbursts were of low am- 
plitude (1—1.5 mag) forcing Kato to classify BF Eri as "a low- 
amplitude dwarf nova". However from the formal side, DNs are 
an subset of CVs, which undergo recurrent outbursts of 2-8 mag. 
The observed "outbursts", including the one observed during the 
2006-Nov set, are of much lower amplitude, which leads us to the 
question: are those events the real DN outbursts rather than some 
kind of flares? 

It is now widely accepted that the reason for a DN outburst 
is a thermal instability in the accretion disc, which switches the 
disc from a low-viscosity to a high-viscosity regime CSmakl [19841 ; 
IOsakill996l ; |Las ota 2001). From the observational side this is (usu- 
ally) accompanied by substantial spectral changes occurring con- 
tinuously during all stages of the outburst. More than half of DNs 
exhibit the transition from emission-line spectrum at quiescence to 
absorption-line spectrum at maximum. For the rest if any emission 
can be seen it is often in the form of emission lines buried in absorp- 
tion troughs jWarneJ[l993:lMorales-Rueda & Marshll2002l) . Quite 
often the spectral line profiles undergo dramatic changes as new 
emission sources arise in the system such as the inner hemisphere 
of the secondary star. This emission is caused by increased irra- 
diation from the accretion regions during the outburst and can be 
seen i n hydrogen and neutral helium lines th rough Doppler tomog- 
raphy iNeustroev. Zharikov. & M ichel 2006). Many systems show 
also strengthening of the He II and C III/N III line emissions during 
an outburst. 

None of these outburst signs was observed in BF Eri. Thus the 
observed flares can hardly be associated with a DN outburst. This 
in turn raises a question on the classification of the system. Ac- 
tually, BF Eri does not also show many other usual observational 
signs of DNs. For example, in spite of the high spectral resolution 
used, all the emission lines show not the double-peaked profiles 
but single-peaked ones. From the above-estimated system parame- 
ters we can expect the double peak separation to be more than 500 
km s _1 , that is more than enough to be resolved. However, though 
the profiles are highly variable, they do not show any signs of the 



double-peaked accretion disc structure. In consequence, a ring-like 
emission distribution, characteristic of a Keplerian accretion disc 
about the WD, is also absent from the Doppler maps. On the other 
hand, the bulk of the Balmer emission is located to the lower left 
quadrant of the tomograms. Both the emission line profiles and the 
appearance of the Doppler maps resemble those in nova-like CVs 
which in turn resemble those typical for old novae at minimum 
light. 

By definition, NLs should not display any DN outbursts. 
Their almost steady brightness is thought to be due to their 
mass transfer rate M exceeding the upper stability limit M cr u. 
These stars are thus thermally and tidally stable. However, 
iHonevcutt. R obertson. & Turned i ll 9981) reported that a significant 
fractio n of NLs displays so-called "stunted" outbursts. IHonevcutt! 
d200lh found many similarities between them and outbursts in or- 
dinary DNs but the former are of much smaller amplitude (0.4-1 
mag). We suppose thus that BF Eri might be a NL system exhibit- 
ing "stunted" outbursts. 

The AAVSO database has almost fourteen hundr ed BF Eri ob- 
servations obtained largely after 1999 (Henden 2007). We have an- 
alyzed the light curve compiled from the 'Visual magnitude' and 
'V magnitude' data and found that the outburst activity of BF Eri is 
quite complex (Fig.QJJ] left). Histogram of the AAVSO data reveals 
at least three types of flashes (Fig. [Toj right). Besides the above- 
mentioned low-amplitude flashes with the peak magnitude ~13.0 
mag and ~ 1 3.5 mag, there were also 8 events with the amplitude 
more than 1.6 mag but always less than 2.3 mag, with V ma x ~ 12. 
Formally, the brightest flashes lie on the bottom edge of the ac- 
cepted range of the outburst amplitudes in ordinary DNs and might 
be such outbursts. The nature of the low-amplitude flashes is still 
unclear; further investigation is needed. 

Finally, we have obtained the periodogram, using the AAVSO 
data (Fig. II It. It shows significant excess power in the interval 60- 
90 d, reflecting a reliable stability of the flashes. 



7 SUMMARY 

We have presented spectroscopic observations of the cataclysmic 
variable BF Eri during its low and bright states. The most important 
results of this study can be summarized as follows: 
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Figure 11. The Lomb-Scargle power spectrum of BF Eri, derived from the 
AAVSO data. 



(i) The orbital period of BF Eri is P orb = 0.270881(3) days. 

(ii) We have shown that BF Eri contains a massive WD (Mi ^ 
1.23Mq). This allows us to consider the system as a SN la progen- 
itor. 

(iii) The spectral type of the secondary star is K3±0.5. 

(iv) We have shown that the secondary in BF Eri is an evolved 
star of the mass of 0.50-0.59 Mq, whose size is about 30 per cent 
larger than a ZAMS star of the same mass. 

(v) BF Eri shows the high 7-velocity (7 = —94 kms -1 ) and 
substantial proper motion. With our estimation of the distance to 
the system (d w 700 ±200 pc), this corresponds to a space velocity 
of ~350 km s~ with respect to the LSR. The cumulative effect of 
repeated nova eruptions with asymmetric envelope ejection might 
explain the high space velocity of the system. 

(vi) We have analyzed the outburst behaviour of BF Eri and 
question the current classification of the system as a dwarf nova. 
We propose that BF Eri might be an old nova exhibiting "stunted" 
outbursts. 
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